Percolating granular superconductors 
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We investigate diamagnetic fluctuations in percolating granular superconductors. Granular super- 
conductors are known to have a rich phase diagram including normal, superconducting and spin 
glass phases. Focusing on the normal-superconducting and the normal-spin glass transition at low 
temperatures, we study he diamagnetic susceptibility x an d the mean square fluctuations of the 
total magnetic moment x' 2 ' of large clusters. Our work is based on a random Josephson network 
model that we analyze with the powerful methods of renormalized field theory. We investigate the 
structural properties of the Feynman diagrams contributing to the renormalization of x an d X^ ■ 
This allows us to determine the critical behavior of x' 1 ' and x^ to arbitrary order in perturbation 
theory. 



PACS numbers: 64.60.Ak, 74.80.Bj, 74.25. Fy 
I. INTRODUCTION 

The discovery of high-temperature superconductivity 
raised hope of an important role of this phenomenon 
in applications requiring high current densities, such as 
power transmission lines and high-field magnets. The 
most severe factor limiting current densities is that all 
practical materials contain defects like impurities, grain 
boundaries and other extended defects. Thus, it is im- 
portant to investigate the role of disorder in supercon- 
ductivity. 

Percolation theory [Q plays a predominant role in the 
study of disordered systems. In the course of the years 
it has been applied to granular superconductors in many 
ways [^| . The diamagnetic properties of disordered com- 
posites of superconducting and non-superconducting ma- 
terials can by studied in terms of a percolation model 
where superconducting grains are located on the sites of 
a (d-dimensional) hypercubic lattice and where Joseph- 
son junctions occupy nearest neighbor bonds with a given 
probability p. In the following we will refer to such a net- 
work as a random Josephson network (RJN). 

The phase diagram of the RJN has a rich structure de- 
pending on the occupation probability p, the temperature 
T, and an external magnetic field B. This phase diagram 
was explored in a seminal work by John and Lubensky 
(JL) |}|. Viewing the unoccupied bonds as normal con- 
ductors, one has a normal phase, a Meissner phase, and 
a spin-glass (SG) phase. For type II superconducting 
materials one has in addition an Abrikosov flux lattice 
phase. For p below the percolation threshold p c there 



are only finite superconducting clusters and hence there 
is no macroscopic superconductivity. For p exceeding p c 
there exists at least one spanning cluster and the phase 
depends on B. As p — > p c at B — one encounters for 
sufficiently low temperatures a transition from the nor- 
mal to the Meissner phase (henceforth, we refer to this 
transition as transition I). This Meissner phase has its 
typical hallmarks, viz., expulsion of magnetic flux and 
a non-vanishing average condensate wave function. For 
B > the system crosses over between the normal and 
the SG phase as p c is approached at low T (this transi- 
tion will be referred to from now on as transition II). The 
SG phase is characterized by a vanishing average conden- 
sate wave function but non- vanishing Edwards- Anderson 
order parameter [I] . 

Early theoretical work on the transition I dates back to 
the beginning 1980's. It was predicted by de Gennes || 
and Alexander || that the configurationally averaged 
diamagnetic susceptibility x diverges as 
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where v is the percolation correlation length exponent 
and t ist the conductivity exponent of the random re- 
sistor network (RRN). A few years later John, Luben- 
sky, and Wang (JLW) H presented a renormalization 
group analysis of the RJN based on a replicated Landau- 
Ginzburg- Wilson (LGW) type continuum model intro- 
duced in Ref. 0. They studied a s well as the mean 
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square fluctuations x °f the total magnetic moment 
for transitions I and II. Their results obtained in a per- 
turbatio n ca lculation to one-loop order support the pre- 
diction (01). Moreover, their calculations resulted in 
^( 2 ) = for transition I and x — for transition II. 
X^ was found to diverge at transition II as 
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where ti is a crossover exponent distinct from t. Roux 
and Hansen Q carried out numerical work to calculate 
ip for d = 2 dimensions. As JLW, Roux and Hansen 
relied on linearized network equations. Their result 
i p = 1.36 ± 0.02 agrees well with the scaling relation 
( 1-lbj ) if the established values for v and t are inserted. 
Wang and Lubensky [^) utilized a low concentration sc- 
ries expansion to determine ip for d = 2. Their result 
ip = 1.21 ± 0.0 3, however, is inconsistent with the scaling 
relation ( 1.1b ). Recently, Knudsen and Hansen (l(J car- 
ried out numerical simulations avoiding the linearizations 
involved in Refs. H), and |§- They obtain \p = 1.2 in 
agreement with the series expansion result of Ref . j|] . 

It appears that a linearization of the network equations 
has a crucial effect in d = 2. To our opinion this is 
plausible because the RJN is intimately related to the 
diluted x-y model [llj]. Hence one should expect that 
vortex excitations become important in d = 2 and that 
a linearized description, corresponding to a spin wave 
approximation, is insufficient in this case. 

Last but not least in our little historical reminiscence 
we quote experimental values. Misra and Misiak mea- 
sured ip = 1.32 @ and <p = 1.45 ||. 

The previous paragraphs indicate that the status of 
diamagnetism in the RJN has several ramifications. On 
this basis it is hard to draw reliable conclusions on the 
nature of the diamagnetism in granular superconductors. 
Further work on this subject seems to be in order. 

In this paper we determine the scaling behavior of x 
and %( 2 ) for transitions I and II by the powerful methods 
of renormalized field theory |Q . Our approach is based 
on the LGW type Hamiltonian introduced by JL. Us- 
ing dimensional regularization and minimal subtraction 
we explore the renormalization of %W and X ■ Upon 
analyzing the general structure of the Feynman diagrams 
contributing to these renormalizations, we derive the crit- 
ical behavior of x an d X^ to arbitrary order in per- 
turbation theory. The outline of this paper is as follows. 
In Sec. [fl] we provide background on the model underly- 
ing our work. At first, we concretize the definition of the 
RJN and sketch its microscopic description. We mention 
the key physical quantities that are implicit in the mi- 
croscopic model, such as the average tunneling current, 
and explain how they might be calculated by using the 
replica formalism. Then we condense the microscopic 
into a mesoscopic model that is represented by a field 
theoretic Hamiltonian Ti. Our final expression for Ti is 
corroborated by a subsequent scaling analysis. In partic- 
ular the irrelevance of a certain coupling associated with 
T 2 is revealed. Then we elaborate on several physical 



quantities that are native in the mesoscopic description. 
We give their definitions and explain how to extract them 
in the replica framework. A brief review of the RJN phase 
diagram concludes Sec. Q. Section [II contains the core 
of our renormalization group analysis. We gather the di- 
agrammatic elements that are the ingredients of our per- 
turbation calculation. Then we take a short detour and 
outline the renormalization and the scaling behavior of 
the order parameter correlation functions. This provides 
some background for our main task, the analysis of x^ 
and x'' 2 ''- Next, we calculate the Feynman diagrams con- 
tributing to x an d X*' 2 ' a t one-loop order. Equipped 
with some intuition about these diagrams we then deter- 
mine their structural properties for arbitrary order in the 
loop expansion. At the end of Sec. Ill we describe the 
renormalization of x^ an d x'" 2 *' i se t U P the correspond- 
ing renormalization group equations and determine the 
scaling behavior. In Sec. IV we give a brief summary and 



concluding remarks. Technical details on the derivation 
of the full Gaussian propagator for the RJN can be found 
in appendix |A|. 



II. THE MODEL 

A. The random Josephson network 

As mentioned briefly in the introduction, a RJN con- 
sists of superconducting grains located at the sites i of a 
rf-dimensional hypercubic lattice. Bonds between nearest 
neighboring sites are randomly occupied with probability 
p by Josephson junctions and, respectively, empty with 
probability 1 — p p| . Each grain is characterized by a 
condensate wave function 



yfp exp(i6* 4 ) 



(2.1) 



Note, that the density of Cooper pairs on each grain is 
assumed to be a constant p and that only the phase Oi 
is allowed to fluctuate. The fixed amplitude approxima- 
tion neglects charging effects due to quantum fluctua- 
tions [lq] . It is justified for grain sizes of the order of or 
smaller than both the bulk superconducting coherence 
length and the London penetration depth for the grains. 

The form of the wave function (2.1) leads within the 
tight binding model to a quantummechanical expectation 
value for the total energy given by 



H = - 



<i,j> 



K it j cos((5ij) 



(2.2) 



Ki t j is a hopping matrix element for the Cooper pairs. 
Here it is a random variable that takes on the value 1 
with probability p and the value with probability 1 — p. 
The sum in Eq. ( |2.2| ) runs over all nearest neighbor pairs 
< i,j >• 



(2.3) 
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describes the phase difference between adjacent sites and can be written as 



A% f j — e 



A- til. 



(2.4) 



Here e* is an abbreviation for 27r/$o with <&o being the 
flux quantum. The line integral is taken over an arbitrary 
differentiable curve from i to j. A is the vector potential. 
We consider the gauge field A being entirely determined 
by a fixed external ma gnetic field and neglect fluctuations 
in A. Expression (2.2) governs the equilibrium statistical 
mechanics of the RJN and represents a Hamiltonian in 
the sense of statistical mechanics. Note the following 
important feature of H: it is invariant under the gauge 
transformation 



0i -*■ 0; + a(i) , 
A A Va, 



(2.5) 
(2.6) 



where a is an arbitrary scalar function of the space coor- 
dinate. 

A fundamental role in the RJN is played by the tun- 
neling currents 



dH 



-Kij sin(6j,j) . 



(2.7) 



Their averages and correlation functions represent inter- 
esting observable quantities. Averaging over thermal de- 
grees of freedom, henceforth indicated by (. . .}t, may be 
discussed via the free energy 



F T = -TlnZ 
with the partition function Z given by 



Z = J P6»exp(-T" 1 i7) 



(2. 



(2.9) 



For convenience we have set the Boltzmann constant 
equal to one. J DO is an abbreviation for J J\ i d9i, where 
the product is taken over all lattice sites. 

Of course, a meaningful characterization of the statisti- 
cal properties of the RJN requires more than just thermal 
averaging. In addition, a quenched average [• ■ -]c over all 
possible configurations C of the diluted network needs to 
be performed. This average can be achieved with help of 
the replica trick, n copies of the network are considered 
simultaneously upon introducing the replicated Hamilto- 



H({5}) = j2ms {a) }) = - Z z E ^cos^j), 

a — 1 a—1 <i.j> 

(2.10) 

where S = (5^, . . . , 5^). With this trick the configura- 
tionally averaged free energy 



F 



-TpnZ] 



(2.11) 



with 



F = lim -F n 

n— >0 n 



F n = -Tin [Z n ] c 



(2.12) 



(2.13) 



The key benefit of this procedure is that the problem of 
averaging In Z is basically replaced by the easier task of 
averaging Z n . From the free energy various quantities of 
interest can be extracted upon taking derivatives. The 
average current is given by 



dF n 



(2.14) 



The derivation of correlation functions in the replica 
framework requires a little caution. Note that 



Cg'fl = - hm d ' Fn 



(2.15) 



splits up into a replica diagonal and a replica independent 
part. Due to the permutation symmetry between the 
replicas C l -"' /3 ' ) is of the form 



r (a,P) _ /nr(l)r , s~i 



1,3 



(2) 



(2.16) 



with 



1 {[(ll j ) T ] c -[{Ii, j ) T ]c}-E c (2.17) 



and 



C%=T 



J {[fc);ic-N T ]c}' ( 2 - 18 ) 



The E c emerging in Eq. ( [2.17 ) stands for the condensa- 
tion energy 



Tie 



(2.19) 



Field theoretic Hamiltonian 



Now we proceed towards a field theoretic model for 
the RJN. The following derivation of a field theoretic 
Hamiltonian is guided by the work of JL. 

Our starting po int i s the observation, that the [Z n ]c 
appearing in Eq. (2.12) can be written as 



[Z* 



V6 exp(-T- 1 J ff eff ) 



with an effective Hamiltonian 



H ef f = - In 



exp 



-T 



MM 



(2.20) 



(2.21) 
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By virtue of the replica approach we can perform the 
quenched average once and for all at this early stage. 
This leads to 



(2.22) 



<«,j> 



with 



K 5 



hi 



1 + wcxp [T^E 



a=l 




(2.23) 



where v — p/(l — p). In Eq. (2.22) we have dropped 
a constant term Nb ln(l — p), where Nb stands for the 
number of bonds in the undiluted lattice. 

Next, we adopt ideas developed by Stephen JlTj in the 
context of the RRN. The idea here is to introduce the 
quantity 



= cxp 



(iX • , 



A^O. 



(2.24) 



As we go along, this quantity will grow into the role of 
an order parameter field. The A appearing in ( [2.24 ) is an 
n-component vector in replica space, A = (A^ , . . . , X- n >). 



The dot product in Eq. (2.24) is defined as A ■ 9 



Y^=\ X^9^ a \ The condition A ^= is imposed in order 
to qualify ip X (i) as an order parameter. That is simply 
because tp^(i) is equal to one and hence, being a constant, 
not capable of sensing a phase transition. The compo- 
nents A^ Q ' are chosen to take on integral values. With 
this choice, the exp(iA • 9) represent a complete set of 
orthonormal functions satisfying the orthonormality and 
completeness relations 

-1- ^ d n 9 exp (it ■ 9) = S U (2.25a) 



and 



L^exp(^.«)^(0) . (2.25b) 



(2 



Based on these relations one can define the replica space 
Fourier transform 

4> (S, i) = ex P (-^ ■ 9j i/3 X (i) (2.26) 

of Note that the condition A ^ on trans- 

forms into the condition 



1 



(2*)" J_ n 



d n 9 



(2.27) 



on (f>(6,i). Thus, <fr(9,i) can be interpreted as a (contin- 
uously indexed) Potts spin fll| . 



Now we expand K(5i t j) in terms of Rewriting 

K{5ij) as 



= j d n SK(sjs(s-S it ^ 
we obtain after a little algebra 
K (4-) = E V'xW^-xO') exp (-it • AiJ) K (] 



(2.28) 



A^O 



(2.29) 



K(X) is defined as the replica space Fourier transform of 
K(S): 



K 



f 



(2tt)" 



d n 5 exp (it ■ 8j K (tj . (2.30) 



In writing Eq. ( [2.29 ) we have dropped a constant term 
K(6). To evaluate Eq. fl2.30j ) further, we insert ( |2.23| ) 
and expand the logarithm. We arrive at 



with 



F 



K 



CXJ 



(-1)' 



v Fi 



1=1 



(2.31) 



27T) 



x exp 



n 

^2 Ux^sM + IT- 1 cos (<5 (q) ) 



(2.32) 



The integral in Eq. ( |2.32j ) can be evaluated in the low 
temperature limit by employing the saddle point method. 
Note that this step amounts to linearization of the net- 
work equations, or in other words, to a spin wave approx- 
imation. For n — > this procedure leads to 



K (X) =t + wX 2 + OIT 2 



(2.33) 



where r = r(p) and w — w(p,T) ~ T are expansion 
coefficients. In Sec. [II C| w e will show explicitly, that 
higher order terms in fl2.33| ) lead to irrelevant contribu- 
tions in the field theoretic formulation, and that, hence, 
their omission is justified. 

Collecting we find the following expression for the ef- 
fective Hamiltonian: 

H <*= E E^cw-xO') 

<ij> A#0 

x exp (-it ■ AiA It + wX 2 + . . . J . (2.34) 

At this stage we carry out a gradient expansion. We have 
to pay regard t o the fact that the effective Hamiltonian 
as given in Eq. (2.34) is invariant under the gauge trans- 
formation 



ip?(i) —> ip?(i)exp iX-a(i) 



A(i) -> A(i)-— V3(<). 



(2.35a) 
(2.35b) 
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Hence, we must keep only those terms in the gradient 
expansion that comply with this invariance. We obtain 



A#0 



x{l + i[6i- (v-ie*\-A(i) 
x |r + wX 2 + . . .} , 



(2.36) 



where the second sum runs over all lattice vectors b^ be- 
tween site i and its nearest neighbors. 

We proceed with the usual coarse graining step and re- 
place V'xW by the order parameter field "0^(x). The or- 
der parameter field inherits the constraint A ^ 0. Then, 
a mesoscopic free ener gy in the spirit of Landau is de- 
vised. Guided by Eq. (2.36) we write down the Landau- 
Ginzburg- Wilson type Hamiltonian 

n = J d d xi I Y, ^-a (*) K ( v > x > a) i>x W 

X#0 

V^ X (x)^_ X , (x)^ +X; (x) 1,(2.37) 



E 

X,X',X+X'#0 



where terms of higher order in the fields have been ne- 
glected since they turn out to be irrelevant. The kernel 
appearing in the Eq. (2.37) is given by 



K 



(v,a,a) 



T + W\ 2 



iX ■ A 



(2.38) 



The coefficients r and w should be understood as the 
coarse grained analogs of the original coefficients featured 



in Eq. (2.36). Similarly, the A in the kernel (2.38) is a 
coarse grained version of the original gauge field. The 
coarse grained A is defined so that it incorporates the 
charge e*. 

It must be emphasized that TL is invariant under the 
gauge transformation 



V^(x) — > ^(x)exp iX ■ a(x) 
A(x) -> A(x) - Va(x) , 



(2.39a) 
(2.39b) 



with the components of a being arbitrary scalar func- 
tions of x. This gauge invariance will have important 
consequences as we go along. 

We point out that TL resembles for vanishing A the 
form of the field theoretic Hamiltonian for the RRN as 
studied by Harris and Lubensky and the present au- 
thors 1 20, For vanishing A, the only formal distinc- 
tion resides in the different domains of A. In the limit 
w —> p2| , however, this difference has no consequence 
and the perturbation expansions for the RJN and the 
RRN coincide. For w = 0, in particular, both models 
reduce to purely geometric percolation. 



C. A note on relevance 

Here we will show that it is indeed justified to truncate 
the expansion (2.33) at first order in T. In other words, 
we will show that the higher order terms are irrelevant in 
the sense of the renormalization group. Our actual tool 
will be a scaling analysis in the replica variable A. 

Now suppose we had retained higher order terms in 
the expansion ( p. 33 ). Then the kernel of TL would be of 
the form 



K 



(v,a,a) 



wX 2 



w k X 

k=2 



2k 



V-iA-A 



(2.40) 



with Wk ~ T ■ We facilitate our scaling analysis by set- 
ting A — » b X, where b is some scaling factor. Upon 
substituting V'xM = ^l-iX^^ ^ e Hamiltonian we 
get 



TL 



^_ lX (x),A(x);r,w,{u; fe } 



^ aVo 



+ - V V f -(x)V f (x) 




(2.41) 



Renaming the scaled replica variables A' = b X A leads to 



n 



^Xt ( x )' A(x);r,w,{w; fe } 



+ 



/ { \ E 4 w K ( v > &xt > A ) (x) 

^ X^O 

X,X',X+X'^0 ' 



(2.42) 



Now we are going to exploit an important feature of the 
summations over the replica variable A. In the low tem- 
perature limit, i.e., for w — * 0, the summation X^X^o ' ■ ■ 
can be replaced by the integration d n X .... Pois- 
son's summation formula guarantees that the neglected 
terms are of the order exp(— const/w) p3fl . In the contin- 
uum formulation the rescaling leads to b n j_ d""X' .... 
Hence, the scaling factor b drops out in the limit n — > 
and we can identify A^ with A. We are led to the conclu- 
sion 



n 



V^-ixW' A(x);r,u),{w fe } 



= TL 



^ x (x),bA(x);T,b 2 w,{b 2k w k } . (2.43) 
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Next we consider the implications of Eq. (2.43) on the 
free energy. In the present field theoretic formulation, 
the Helmholts free energy is defined as 



T 



A(x);T,T,w,{w k } =-T\nZ (2.44) 



with the partition function 



Z = fv^ exp (-T^H [^ x (x), A(x);r, w, {w k } 



(2.45) 

Here, T>i\) indicates an integration over the set of variables 
{■0x( x )} f° r au x an d A. Equation ( 2.43| ) implies that 



A(x);T, r,w,{w k } 



= T 



bA(x):T,T,b 2 w,{b 2k w k } 



(2.46) 



Of course, we are free to choose the scaling parameter to 
our liking. With the choice b 2 = w^ 1 we obtain 



T A(x);T, r, w, {w k } 



= T 



■^WiT,r,l,{3}]. (2.47) 



We learn from Eq. ( 2.47 ) that the coupling constants 
w k appear only in the combination w k /w k . A trivial 
consequence of the fact that the Hamiltonian Ti must be 
dimensionless is that wX 2 ~ fi 2 and w k X 2k ~ fi 2 , where 
/i is an inverse length scale. In other words, w\ 2 and 
w k X 2k have a naive dimension 2. Thus, Wk/w k ~ /i 2 ~ 2fe 
and hence the w k /w k have a negative naive dimension. 
This leads to the conclusion that the w k /w k are irrelevant 
couplings and that the leading critical behavior of the free 
energy is described by 



T [a(x); T, t,w\ = / [w~ 1/2 A(x); T, r 



(2.48) 



where / is some functional of u; _1 / 2 A(x). 

The fact that W2 appears only in an irrelevant com- 
bination was overlooked in Ref. |Q. This ultimatively 
led to an erroneous prediction for the scaling behavior of 

v( 2 ) 



D. Current density, magnetization and related 
quantities 

In this section we elaborate on various physical quan- 
tities embedded in the field theoretic model. We provide, 
within the replica framework, definitions of the current 
density and the magnetization along with their averages 
and correlation functions. We explain the physical con- 
tent of replica quantities and describe how it can be ex- 
tracted. 



The role occupied in the original microscopic model by 
the replicated tunneling currents I^ a ) is taken in the field 
theoretic formulation by the replicated current density 

Joe), \ SH 



(2.49) 



Note that the index i specifies here the component of the 
current density in d-dimensional space and should not be 
confused with the site i. 

The current density has a very important feature. It 
represents the Noether current associated with the gauge 
invariance of the Hamiltonian Ti. Hence, it satisfies the 
conservation relation 



V- J (Q) (x) = 



(2.50) 



Averages in the field theoretic formulation are declared 
by means of the functional integral 



= 1 J ^•••exp(-T- 1 W) 



(2.51) 



Various of these averages can b e extracted from the free 
energy introduced in Sec. II C. This free energy can be 
expanded as 



T A =^ - 



d d x (j( Q) (x))A (Q) (x) 
\J d d x J dVGf^Cx-xOA^W^fCx 



+ ..., (2.52) 

where the summation convention is understood for the 
indices labeling space and replica coordinates. 



t\*)) = - 



5T 



K (Q) (x) 



(2.53) 



A=0 



is the average replica current density. The second order 
term features the correlation function 



C 



(x-x') = 



5 2 F 



A=0 



= T 



1 {(4 a) (x)jf(xo)-(4 a) (x))(jf(xo)} 

- 5(x - x')<^ Yl Xia)xW (V'aWV^x')) • (2.54) 

The relation between the field theoretic average of the 
replica current density and the thermal and configura- 
tional average of the physical current density is straight- 
forward: 



hmU a) (x) 

n— >0 



[(^(x))] 



c 



(2.55) 



As far as correlation functions are concerned, the sit- 
uation is somewhat more subtle. The structure of the 
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correlation functions must be so that they are invariant 
under permutations of the replicas. Hence the two-point 
functions are of the form 



Ct 0) (x - x') = (x - x') 6 a , + Cg (x - x') 



1,3 



(2.56) 



with 



Cl ! y(x-x')=T- 1 {[(J l (x)J J (x')) T ] c 
- [<J<(x)) T (J,(x')) T ] c } - <M(x - t'JBW 



(2.57) 



and 



Cg(x-x') =T- 1 { [(J 4 (x)) T (J,(x')) T ] c 

- [(Ji(x)) T ] c [( J,(x')) T ] c } - <M(x - x')^ 2) , 

(2.58) 

where the replica limit n — > is understood. .Ec and 

(2) 

are the replica diagonal and the replica independent 
part of 

= £ A (-) A (« ^(x)^(x')} • (2.59) 

A>0 



From Eq. (2.57) one learns that 



C«(x-x') = [C iJ (x,x') 



(2.60) 



for 



n — > 0, i.e., the replica diagonal part of c}" (x,x') 
corresponds to the average of the physical density current 
correlation functions Cij(x,x') over all configurations C. 

(2) 

The physical content of C\! will become clear below. 

Now we shift focus and turn from the current densities 
to the magnetization and its correlations. An external 
magnetic field can be introduced into the model via 



F itj (x) = d i A j {x)-d j A i (x). 



(2.61) 



For d — 3 this reduces to the usual B = rotA with the 
components of the magnetic field given by B% — F1.2 and 
so on. Having Fi t j at hand, we can rewrite the expansion 
of the free energy as 

= F[3\ - J d d x (J| a) (x))f(x) 

~\Jd d *J rfVcW) (x _ x ^W( x )^) (x ') 
+ .... (2.62) 
In recasting the second order term we have exploited that 

(2.63) 



by vi rtue of the ga uge i nvariance and its manifestation 
(2.50). Due to Eq. (2.63), the Fourier transform 



Cg'^k) = J d d x Cg ,/3) (x) exp (-ik • x 



) (2-64) 



of cf^'^Vx) is of the form 



C 



1, j 



(k) = (k 2 ^ - hkj) (k) . (2.65) 



The C<- a >fi(x) in the expansion ( [2.62| ) is defined as noth- 
ing but the Fourier transform of C^ a '^(k). 

In the remainder of this paper we will be concerned 
with a homogeneous external magnetic field perpendic- 
ular to the 1-2-plane (x-y-plane). In three dimensions 
this corresponds to a B field pointing in the 3-direction 
(z-direction). We formulate the envisaged magnetic field 
by setting 



with 



Fk,% = B (4,1^,2 ^ 4,2<So) 



(2.66) 



(2.67) 



With help of this expression we obtain the free energy as 
a function of the replicated field amplitude B: 

T(B) = F(6) - J d d x M< Q ) (x) 

~\Jd d xJ A'x W) (x-x')Bt°)BM 



(2.68) 



where 



M<«>(x) = - [( J 2 (Q) (x)) xx (Jr(x)) x 2 \ (2.69) 



and 



X (Q,/5) (x-x') = C (a '^(x-x'). 



(2.70) 



Adapting the usual definition of the magnetization in 
the homogeneous field setup to the replica framework we 
have 



V dB^ 



B=a 



( 2 - 71 ) 



where V stands for the volume of the system and 
is the replica version of the total magnetic moment 

Mtot - J d d xM( X ). (2.72) 

The physical magnet izati on M is retrieved by taking the 
limit n — > 0. From (2.6E) one obtains immediately that 
M vanishes for B = 0. Turning to the diamagnetic sus- 
ceptibility we have 



>,/3) _ 



d 2 T 



1 



V dBWdBW 



B=0 



= i j d d x J d d x' X {a ' P) ■ 



(2.73) 
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A glance at ( 2.7C| ) brings about two important observa- 
tions. First, we see that 



(2.74) 



This relation will play an important role in our actual 
calculations. The second observation is that the diamag- 
netic susceptibility should have the same replica struc- 
ture as the current density correlations, i.e., 

X K/3) =X (1) ^ + X (2) . (2.75) 
In the replica limit x has the physical content 



X (1) = [X\c 



(2.76) 



with x being the diamagnetic susceptibility for a given 
configuration C. x^ contains the fluctuations of the 
total magnetization. For n — > one has 



X 



(2) 



f 

V 



d a x I d d x' 



x T-i{ [M(x)M(x')] c - [M(x)] c [M(x')] c } 



= T- 1 {[Ml t ] c -[M tot } 2 c } 



(2.77) 



In the following we will refer to x^ an d X^ m a brief 
fashion as susceptibilities. 




FIG. 1: Schematic phase diagram of the RJN. The critical 
surface (1) separates the insulating and the superconducting 
phase. The critical surface (2) lies between the insulating 
and the SG phase. Both surfaces merge in a line of bicritical 
points (solid) where the insulating, the superconducting, and 
the SG phase meet. For w — there is a critical line (dashed) 
separating the insulating and the SG phase. At the origin 
there is a critical point (dot) between the three phases. 



E. Review of the phase diagram 

In favor of a self contained presentation, we now briefly 
review the phase diagram of the RJN Jjl , see Fig [|. In 
mean-held theory, the phase diagram can be mapped 
out by determining those combinations of the parame- 
ters p, T, and B for which the Gaussian part of Ji devel- 
ops the eigenvalue zero. For convenience, we write the 
space coordinates x = (a;,y, xj_), where xi lies in the 
(d — 2)-dimensional subspace perpendicular to the x-y- 
plane. In the following we use the Landau gauge, i.e., we 
set A(x) = (0, x, 0j_). The eigenvalues of the Gaussian 
part (the Landau levels) can be determined by standard 
textbook methods. One finds 

E(q ± ,X,m) =r + q 2 1 +ioA 2 + (2m + l)|w(A) | (2.78) 

The momentum is the Fourier transform of xj_. m 
labels the Landau levels and takes on the values m = 
0, 1,2,.. .. u>(X) is a cyclotron frequency given by 

u(X)=B-X. (2.79) 

In the following we will assume that the external mag- 
netic field is replica symmetric and write B = for 
every a. 

Obviously, the lowest eigenvalue is associated with 
q ± = 0j_. Now suppose that w > and B > 0. By in- 
spection one finds two modes associated with vanishing 
eigenvalues: (i) One with Ai = (1, 0, . . . , 0) and m = 0. It 



is soft, i.e., its eigenvalue vanishes, for t + w + B = 0. (ii) 
One with A2 = (1, — 1, . . . , 0) and arbitrary m that is soft 
for t + 2w = (note, that t + w + B = t + 2w for B = w). 
Ai and X 2 should be understood as representative for a 
set of equivalent (i.e., identical up to a permutation of 
the components) replica vectors. We identify two crit- 
ical surfaces. The surface specified by r = — (w + B) 
and B < w separates the insulating and the supercon- 
ducting phase. The order parameter for this transition 

is (tpx ( x )) [(^ 6ic )t]c- The surface corresponding to 
t = — 2w and B > w separates the insulating and the 
SG phase. Here, the order parameter is of the Edwards- 
Anderson type, (^ 2 (x)) ^ [( e ^) T ( e -^) T ] c . The 
two critical surfaces merge at a line of bicritical points, 
given by t — ~2w = —2B, where the insulating, the 
superconducting, and the SG phase meet. 

Now we come to the transitions I and II that are the 
main concern of this paper. Admitting w — 0, one finds 
a critical line specified by r = w = 0. Crossing over 
this line by tuning t > (p < p c ) to r < (p > p c ) 
one gets from the insulating to the SG phase. This is 
our transition II. The order parameter here is (^>?(x)) 

with all A =/= satisfying 53«=i = 0- Finally, there 
is the critical point r = w = B = that represents the 
terminus of the two critical lines. Tuning t > to t < 
about this point one crosses over from the insulating to 
the superconducting phase. This is our transition I. Its 
order parameter is (^/jj(x)) with arbitrary A 7^ 0. 
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III. RENORMALIZATION GROUP ANALYSIS 

A. Diagrammatic elements 

To set up a diagrammatic perturbation expansion we 
need to identify the elements contributing to our Feyn- 
man diagrams. Evidently there is the vertex —g/T. 
The Gaussian propagator for the present problem is not 
straightforward to determine, at least not in a non- 
approximate and closed form. That is due to the pres- 
ence of the gauge field. A possible approach is to expand 
the propagator in terms of the c igenfunctions belonging 
to the Landau eigenvalues ( 2.78 ). This route was taken 
by Lawrie |24j in studying the related problem of the 
LGW superconductor. By summing over all Landau lev- 
els Lawrie brought the LGW propagator into an elegant 
form that made multi-loop calculations tractable. In Ap- 
pendix [A] we present a simpler approach that allows us 
to derive the propagator directly without resorting to a 
expansion in terms of Landau levels. Yet, our approach 
reproduces the closed form found by Lawrie. We obtain 



G bold (x,x',A,) =Texp 



x J G(k,A)exp[ik- (x-x')] (l - <%,ff) (3.1a) 

as the principal propagator for the RJN. Here, J k is the 
usual shorthand notation for l/(2-K) d J d d k. G(k, A) is 
given by 



G(k, A) 



x exp 



ds 



o cosh 



(r(A) 



sh (w(A)s) 

tanh ^cj(A)s^ 



(p 2 + q 2 ) 



(3.1b) 



where p and q are conjugate to x and y, respectively. t(A) 
is a shorthand notation for r + wX 2 . The factor (1 — 5^ ~) 

implements the constraint A 0. 

We annotate that one could discuss the RJN phase 
diagram by analyzing the infrared behavior of the prop- 
agator^ 3.1) instead of the minima of the Landau lev- 
els ( 2.78 ). Basically, one just has to consider the limit 
s — » oo and to determine the parameter combinations 
for which the propagator becomes long range. Of course, 
one finds the phase diagram discussed in Sec. II E 



We observe that the principal propagator decomposes 
into two parts, 



G bold (x,x',Aj = G cond (x, x' , A) - G ms (x - x' 



(3.2) 



One of them 

G cond (x,x',A) =Texp 



i^-( X + X ')( y -y') 



< J G(k, A) exp [ik ■ (x - x')] , 



(3.3) 



carries A. The other one, 



G" 



= T 



S x,d 



k T 



exp [ik ■ (x - x')] , (3.4) 



does not. The notation we use here reflects the close 
analogy to the RRN, where A plays essentially the role of 
a current. For the RRN the decomposition of the prin- 
cipal propagator culminates in a real world interpreta- 
tion fl2C " 



\2% 25 



26f £7j £8|, |2JJ £0|, £1|, |3J, |33j £4j, |3j 

in that the Feynman diagrams are viewed as conducting 
networks composed of insulating and con ducting propa- 
gators. The insulating propagator (3.4) is identical in 
form to its counterpart for the RRN. The conducting 
propagator (3.3) reduces formally to its analog for the 



RRN for vanishing cyclotron frequency, 



G 



bold/ 



C,X, A) 



T I exp [ik ■ (x - x') 



k t + w\ 2 + k 2 



= G(x-x',A) 



(3.5) 



Comparing G(x — x', A) to the full conducting propaga- 
tor G cond (x, x', A) it is apparent that the perturbation 
theory simplifies tremendously for vanishing lu(X). This 
simplification will allow us to study the susceptibilities 
at transition I and II with reasonable effort. 



B. Order parameter correlation functions 

Here we will discuss the renormalization and the scal- 
ing behavior of the order parameter correlation functions 

G N ({x, A}; t, w, B, g) = ( Xl ), . . . , ^ N (3.6) 

where we drop the redundant scaling variable T for nota- 
tional simplicity. Though these correlation functions are 
not the main concern of this paper, they deserve some 
attention. First, they are interesting in their own right. 
Second, their discussion will provide some background 
for the subsequent analysis of the susceptibilities. Most 
of the techniques we are going to use, such as dimen- 
sional regularization and minimal subtraction, belong to 
the standard repertoire of renormalized field theory, cf. 
Rev. || . 

To remove ultraviolet (UV) divergences from the order 
parameter correlation functions finite, we use the renor- 
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malization scheme 

ip^tp = Z 1/2 iP, (3.7a) 

r->f = Z- 1 Z t t, (3.7b) 

w — > u> = Z^Z^ui , (3.7c) 

- Z-^Zy^,- 1 / 2 ^/ 2 ^ 1 / 2 ^ 72 , (3-7d) 

B —> B = Z- 1/2 Z B B (3.7e) 

where the ° indicates unrenormalized quantities, /i is the 
usual inverse length scale, e = 6 — d specifies the de- 
viation from the upper critical dimension 6. The factor 
G e = (47r)- rf / 2 r(l + e/2) is introduced for later conve- 
nience. 

At first, we consider the role of the magnetic field. For 
the closely related problem of the LGW superconductor 
it was demonstrated by Lawrie |^| in a two-loop calcu- 
lation that the B-field does not require renormalization, 
i.e., that B = B up to second order in e-expansion, and 
that the remaining Z- factors are independent of B. One 
can show, however, that the validity of these findings is 
not limited to second order; they are valid to arbitrary or- 
der in e-expansion. We will address these points in some 
detail in a forthcoming publication ]37j on the LGW su- 
perconductor. The quintessential points can be sketched 
as follows: the non-renormalization of B is a consequence 
of the gauge invariance. Exploiting this invariance, one 
can show that the current density renormalizes trivially. 
Because the vector potential is conjugate to the current 
density, it does not require an independent renormaliza- 
tion factor also. In turn, B renormalizes trivially, i.e., 
Z~ x / 2 Zb = 1 to arbitrary order in e-expansion. The B- 
independence of the Z-factors follows from the fact that 
B is not dimensionless at the upper critical dimension. 

For the present problem the reasoning of Ref. [[[tJ has 
to be modified somewhat. That is because the magnetic 
field appears in case of the RJN always in the combi- 
nation w~ x l 2 B ~ fi as opposed to the pure B ~ [i 2 
in case of the LGW superconductor. As a consequence, 
the mass r requires in the present problem an additive 
renormalization proportional to B 2 . This subtlety has 
no consequence for our main results and we will ignore it 
in the following. 

The most economic way to determine Z, Z T , Z w , and 
Z g is to exploit the close relation of the RJN to the RRN. 
As indicated earlier, the corresponding two diagrammatic 
expansions coincide for B = in the replica limit. Fur- 
thermore, the RRN reduces for w — to purely geomet- 
rical percolation. Hence, Z, Z T , and Z g are nothing but 
the usual percolation Z-factors known to third order in 
e-expansion pS| ). Z w may be gleaned to second order in 
e from our work on the RRN po|, ^l| . 

Having determined the Z factors, we are now in the 
position to infer the scaling behavior of the order pa- 
rameter correlation functions from their renormalization 
group equation (RGE). This RGE is a manifestation of 
the fact that the unrenormalized theory has to be inde- 
pendent of the arbitrary length scale introduced by 



renormalization. Hence, the unrenormalized correlation 
functions satisfy the identity 



-7^-G N ( jx, A j ;i,w,B,g) = . 



Equation (3.8) translates via the Wilson functions 



7... (u) 

K (u) 

C(u) 



9 1 7 

Mtt - ln z .. 
8/i 

8u 



8 ft 



o 

91nr 



8fi 

8 law 



u (37 - 7„ - e) 

= 7 - 7r , 

= 7 - 7w , 



(3.8) 

(3.9a) 
(3.9b) 
(3.9c) 
(3.9d) 



(the |o indicates that bare quantities are kept fix while 
taking the derivatives) into the RGE 



8 n 8 8 8 N 

8/i 8u 8t 8w 2 



(3.10) 



The RGE can be solved in terms of a single flow pa- 
rameter £ by using the characteristics 



M 
di 




= M , 




(3.11a) 


8u 
di 


= (3(u(i)) , 


«(1) 


= w, 


(3.11b) 




= «(«(*)) , 


f(l) 


= r, 


(3.11c) 


8 

i-\n W 


= C(fi(0) . 


w(l) 


= w , 


(3. lid) 




= l{u{i)) , 


2(1) 


= 1. 


(3. lie) 



These characteristics describe how the parameters trans- 
form if we change the momentum scale /i according to 
fi — > fi(i) = ifi. Being interested in the infrared (IR) 
behavior of the theory , we study the limit £ — > 0. Ac- 
cording to Eq. ( 3.11b| ) we expect that in this IR limit 
the coupling constant u(i) flows to a stable fixed point 
u* satisfying (3(u*) = 0. The IR stable fixed point so- 
lution to the RGE is readily found. In conjunction with 
dimensional analysis (to account for naive dimensions) it 
gives 

x G N (fix, X} ; £- l / v T , e-* /v w, r 2 B, u*,^ 

(3.12) 

with the critical exponents for percolation rj = j{u*) and 
v = [2— k(u*)] _1 known to third order in e |3q |. cf> = v[2— 
C{u*)] is the percolation resistance exponent known to 
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second order in e [£0|, ElJ f39| . We have not yet exploited 
the freedom to choose t. By setting for example I = \t\ u 
we find that the order parameter correlation functions 
scale like 

G N ({x,\};t,w,B,u,») =C {d - 2+vW2 

x n N ({r x x,^ /2 ^ 1/2 A} ; e 2 - 0/2 ^- 1/2 i?) , 

(3.13) 

where £ is the correlation length and Qn is a scaling 
function. 



C. Susceptibilities 

In this section we discuss the renormalization of the 
susceptibilities as well as their scaling behavior for tran- 
sition I and II. We start by analyzing the Feynman di- 
agrams contributing to the susceptibilities. First, we 
consider the one-loop order in some detail. More im- 
portant than yielding concrete results, this provides us 
with some intuition about the general structure of the 
diagrams. Then, this general structure is analyzed for 
diagrams with an arbitrary number of loops. We demon- 
strate how to renormalize the susceptibilities properly. 
Finally, we derive their scaling behavior. 



1. Diagrammatics: one-loop calculation 

Our starting point here is the definition of the current 
density correlation function c[^'{x — x') in Eq. ( p. 54 ). 
At first, we express this correlation function directly in 



terms of the order parameter field. Recall that C\ ■ (x- 
x') is defined originally in terms of the current density 

4 a) (*) = £4 A(Q) 



A^O 



x (x) diil> x (x) - ^ (x) c^_ x (x) ] , (3.14) 

where we have set A = 0. Being a composite field, the 
current density is inconvenient to handle in actual cal- 



culations. Hence, we substitute (3.14) into the current 



density correlation function. Then C|"'^(x — x') is com 



posed directly of correlation functions of the order pa- 
rameter field and its derivatives. Next, we expand the 
weight exp(— T _1 7i) in powers of the coupling constant 
g/T. At zeroth order in g/T we obtain 



a 



(«,/3) 



j T -i [G-boid^ x X) ^G bold (x, x ', A) 

d t G hold (x, x', A) d'j G bold (x' , x, A) ] 
<5(x-x')^G bold (0,A)}. 



(3.15) 



Here it is understood that the bold propagator is evalu- 
ated at w(X) = since we focus on transition I and II. 
Now we decompose the bold propagator into its conduct- 
ing and its insulating part. We recall that the insulating 
propagator contains a factor 5 X g. Due to the A^A^' 
in ( |3.15| ) the all terms containing insulating propagators 
drop out. We obtain upon sending x' — * 



^i_ A («) A («{ T -i[ G ( X5 X)^G(x,A) 

A 

-d l G(x,X)d J G(x,X)] -<S(x)(5 4J -G(0,A)}. (3.16) 

Next we switch to momentum space. Applying the 
usual Fourier transformation to Eq. ( [3.16 ) yields 



A ,/q 
J 2 QiQj 



[r(A)+(q + k/2)2] [r(A) + (q - k/2)>] 

5, 



J i,3 



(3.17) 



Since we are ultimately interested in the susceptibilities, 
we should look at & a ^\\n) rather than G^ ,/3) (k). Via 



taking the trace on both sides of Eq. ( 2.6 5| ) we find that 
2T 



CT(a.ffl(k) = 



[r(A) + (q + k/2)2] [r(A)+(q-k/2)2] 
1 ] 



[r(A)+q^ 2 j 

Here we have used that 
d 

"(A)+q 2 



q [^W+q 2 ]' 



(3.18) 



(3.19) 



in di mens ional regularization. By virtue of the rela- 
tion ( [2.74 ) we obtain upon expansion in the external mo- 
mentum k 



T 



[-(A)+q 2 ] : 



(3.20) 



This is the replica susceptibility at one-loop order. 

We proceed with simplifying the summation over A. In 
Schwinger representation, this summation is of the form 



exp — sr(A) 



(3.21) 
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At this point, we find it convenient to introduce the short- 
hand notation 



I \ = 



. exp 



"(A) 



Revisiting Eq. ( p. 75 ), we deduce that 

(\^\W)x = a5 a , l3 + b, 



(3.22) 



(3.23) 



where a and b are coefficients that need to be determined. 
By analyzing the cases a = /3 and a ^ f3 we find 



1 



b = 



n(n — 1) 
1 



n(X 



n(n — 1) 



n 

E^ (Q 



(3.24) 
(3.25) 



To evaluate a and 6 further we now look at transition I 
and II separately. At the transition I we have B = 0. 
For B = 0, the system is, like the RRN, rotationally 

invariant in replica space and hence ((X^"=i A*-"- 1 ) )a = 
(A 2 )a- Consequently, we have a — -^{X 2 )\ and 6 = 0. 
This leads for the susceptibilities at transition I to 



E 



q KA)+q 2 ]' 



X (2) =0. (3.26) 



At transition II we have 2a=i = ^ anc ^ thus a = 
O(n ) and b = ^(A 2 )a- For the susceptibilities at transi- 
tion II this leads to 



0, X 



(2) _ 



T 
6n 



E 



A 2 * 



q [^(A)+q 2 ] 2 



(3.27) 



Now we are in the position to carry out the remaining 
summations over A along with the momentum integra- 
tion. We outline the remaining steps at the instance 
of a t transition II. Implementing the constraint 
\<.«) = via the integral £ d9 e l9 £" =1 we 

have 



/(2) 



T 
6n 



dee 



q [r(A)+q 2 ]- 



(3.28) 



We find it convenient to use the Schwinger representation 
and rewrite x as 



/(2) 



x exp 



T 
6n 



ds 



xwX 2 — s 



2tt _ 
+ w\ 2 



dB 



E 



q 2 ) +i8j2* {a) 



z=0 

(3.29) 



At this point it is useful to exploit Poisson's summation 
formula and approximate the summation by the in- 
tegration d n X. This integration is Gaussian and it 
can be evaluated straightforwardly by completing squares 
in the exponential. Next, we carry out the particularly 
simple momentum integration. Since we eventually have 
to take the replica limit n — > 0, we then expand x^ m 
powers of n. Upon taking the derivative with respect to 
z and carrying out the integration over 9 we get in the 
replica limit 



X 



(2) 



T 
12 



ds ■ 



1 



6sw 2 



(3.30) 



( 47rs )d/2 

The remaining integration over the Schwinger parameter 
s represents no difficulty. For the renormalization group 
treatment that we have in mind we will only need the 
UV divergent parts of x- These are extracted in form 
of e-poles by expanding \ m powers of e. This provides 
us finally with the following divergent parts of for 
transition II: 

Tr 2 G e f 7r 2 t 

-< 1 -I 

w 12e| 18 w 



(2) 



(3.31) 



Note that the 1-loop result for \dil gi ven m Rev. M is 
erroneous. It incorrectly features W2- 

for transition I can be calculated in an analogous 
manner. We merely need to replace the integration that 
enforces the constraint $^q=i X( a > = by unity. This 
leads to the result 



-Adiv 



Tt 2 



12 e 



(3.32) 



At this point a comment on the r dependence of the sus- 
ceptibilities is in order. From a technical point of view, 
dimensional regularization is the most convenient way of 
dealing with UV divergences. However, dimensional reg- 
ularization has some unphysical features which are inti- 
mately related to its simplicity. In the less economic but 
more physical cutoff regularization one treats UV singu- 
larities by introducing a cutoff A and by replacing the full 
integration J q by /| q |<^- This procedure leads typically 
to terms proportional to In A and terms varying as some 
power of A. The logarithmic divergences for A — > oo 
have their analog in dimensional regularization in form 
of the £ poles. The terms algebraic in A, however, are 
unaccounted for in dimensional regularization. In case of 
the susceptibilities this neglect conceals essential physics. 
Hence, the terms algebraic in A must be incorporated. 
We get 



(i) = T 

Xdiv w 



A^ + A^r-^-r 2 
12 e 



(3.33) 



and 



,(2) _ T 



Adiv 



A^ + A^t-^-t 2 

12 £ 



T 



Bo A + £?i A r + Bitx t 



2 JZ 



7T 2 G, 



5 r 3 



12- 18s 



(3.34) 
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FIG. 2: General structure of the Feynman diagrams con- 
tributing to the susceptibilities. The hatched blob symbolizes 
an arbitrary number of closed loops composed of the vertex 
—g and conducting and insulating propagators. To the left 
and the right of the blob we have conducting propagators. 
The wavy lines indicate insertion points for external momenta 
into these propagators. The two bars stand for a factor A ■ A'. 



where the ^4's and B's are numerical coefficients. 



2. Diagrammatics: general structure 

Higher orders in the expansion of the the weight 
exp(— H) in powers of the coupling constant g correspond 
to multi-loop diagrams. In Schwinger representation, the 
A-featuring part of such a diagram is of the form 



A («) A '(/3) 



exp 



p 



(3.35) 



Here, {k} stands for a complete set of independent loop 
currents. The sums over p are taken over all conduct- 
ing propagators. The currents A p running through the 
conducting propagators are linear functions of the loop 
currents, X p — X p ({k}). The currents A and A' are iden- 
tical to particular A p . In the one- loop example given in 
the preceding section, there is one loop current A and the 
two A p as well as A' are identical to A. 



Of course, ( 3.35 ) splits up into a replica diagonal and 
a replica in depende nt part. Generalizing the argu ment s 
given in Sec. Ill C 1 it is not difficult to deduce from ( 3.35| ) 
that vanishes at transition I to arbitrary order in 
the loop expansion. The same goes for x & t transition 
II. Moreover, one finds that the A-featuring part of non- 
vanishing diagrams is of the structure 



— A • A' exp 

{£} 

1 d 



nw dz 



E 



exp 



£ v^-^]TaW 

p \ a—l 

zwX ■ A' 



(3.36) 



2=0 



In case of the transition I the P are all zero. Figure || 
depicts these non- vanishing diagrams. 



Now to the evaluation of Eq. (3.36). Its right hand 
side factorizes into n equivalent factors. Exploiting the 
fact that the X p are linear functions of the loop currents, 
we write 



1 d 



nwdz 



n e 

a=l = 

*5>({0})k 



exp 



LI' 



(a) 



(3.37) 



2=0 



for the right hand side of Eq. ( |3.36| ). The sums over I and 
I' run over the complete set of independent conducting 
loops corresponding to {/?}. Aiji({s} 7 z) and bi({8}) are 
linear functions of their variables. In case of the tran- 
sition I the h({9}) are all zero. Since we are interested 
in the limit w — > 0, we may apply Poisson's summa- 
tion formula, i.e., replace the summations by 

integrations f^^dn^ 01 '. The so obtained integrations are 
Gaussian and hence straightforward. They yield 



1 d 



' dz 



det(wA) 



: exp 



■— b T A-H 

4w~ = 



.(3.38) 



2 = 



Here, b stands for a column matrix constituted by the 
bi {{0}) and A stands for a square matrix with the ele- 
ments Ai ji (fs j, z). We extract the limit n — > by ex- 
panding ( |3.38 ) in powers of n and find, up to a minus 
sign, 



d 



2w dz 



-^-k T A- 1 b + In det(io A) 



2w~ 



- 1N , i (det(A))' 

= Afe U - 1 ) k+^~ —j ji! ■ (3-39) 
Aw 2 y = ' 2w det(A) 



In writing Eq. ( 3.39 ) we use the shorthand notations 



oz— 
d 



2 = 



( det G4))' = ^ de1 



(3.40) 



The first term on the right hand side of Eq. (3.39) is 
a homogeneous function of the Schwinger parameters s p 
of degree —1, i.e., is behaves generically like s~ 1 . The 
second term is homogeneous in the s p of degree 0, i.e., it 
goes generically like s°. 

Now to the momentum-featuring part of our typical 
diagram depicted in Fig. ^|. Suppose the diagram has L 
closed loops. Then the momentum integrations result in 
a factor that is a homogeneous function of the Schwinger 
parameters of degree —Ld/2. In other words, the mo- 
mentum integrations yield a generic factor s~ Ld / 2 . 

Next we turn to the integration over the Schwinger 
parameters. Suppose our diagram has P propaga- 
tors. Having each propagator represented in Schwinger 
parametrization, we have P integrations over Schwinger 
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parameters. These may be viewed as a generic factor 
s p . Collecting, we find that the term proportional to 
in Eq. ( 3.39| ) goes like S p - Ld / 2 and that the term 
proportional to w~ 2 goes like S p - Ld l 2 -' V , By inspection 
one finds that the following topologic relations apply to 
our typical diagram: L = P — V — 1 and 3V + 4 = 2P, 
where V denotes the number of vertices. For the w^ 1 - 
term these relation lead to P — i| = — ^=^V + 4— d, i.e., 
at d = 6 this term goes generically like s~ 2 . Similarly, 
one finds a s~ 3 -behavior for the w _2 -term. Now we have 
sufficient information to single out the r-dependence of 
the terms. By a change of variables of the type s — > s/r 
we learn that the w _1 -term is associated with a factor t 2 
whereas the w _2 -term features a r 3 . 

For transition II we still have to deal with the inte- 
grations over the (9's. These, however, merely result in 
purely numeric factors. Harvesting the findings of the 
above reasoning, we deduce that the divergent parts sus- 
ceptibilities have the general structure 

■X(u), 



-^div 



(2) 



W 

Tt 3 



Tt 2 

, F ti j X(u). 

z w 



(3.41) 
(3.42) 



It is understood that Eq. ( 3.41 ) and Eq. ( 3.42 ) refer to 
transition I and transition II, respectively. The coeffi- 
cients X(u) and Y(u) have a Laurent expansion of the 
form 



X(u) 



X k (u) 



fc=i 



with 



X k (u)= J2 Xk 



i>k-l 



Y k (u) 



with Yfe (u) 



fc=i 



E 

i>k-i 



ft 



(3.43) 



(3.44) 



with Xk, m a nd ft. m being nume rical c oeffic ients. As ar- 
gued in Sec. |III C l| , Eqs. ( |3~4l"l ) and (|1|) have to be 
supplemented by terms varying as powers of A. This 
gives finally 



y (1) - 

Adiv 

(2) . 
Xdiv 

T 



T 

w 
T 



A,(ii)A 4 + A^A 2 t + X(u) t 2 
A (u)A i + yli(ii)A 2 r + X(u) r 2 



B {u)A 6 + Bi (u)A 4 r + B 2 {u)A 2 



(3.45) 



Y(u)t 3 
(3.46) 



where Aq(u), Ai(u) and so on are u dependent coeffi- 
cients. 



3. Renormalization and scaling 



In Sec. [II B we argued that the order parameter corre- 



lation functions can be renormalized by the renormaliza- 



tion scheme (3.7). This scheme, however, is not sufficient 



to renormalize the su sceptib iliti es. To r emove the e-poles 



encountered in Sees. Ill C 1 and III C 2, one has to resort 



to an additive renormalization 
1 



Tir, 



TLr. 



d d xF id {x)F id {x). (3.47) 



Here 7i re n stands for the renormalized Hamiltonian ob- 
tained form the original bare Hamiltonian Ti. by apply- 
i ng th e renormalization scheme (3.7). We point out that 
( |3.47 ) is not adequate to renormalize the free energy com- 
pletely. As power counting shows, one also needs additive 
counter-terms of third and fourth order in Fjj as well as 
counter-terms containing derivatives of Fij . However, 
for our central task, i.e., for determining the scaling be- 
havior of the susceptibilities, it is sufficient to consider 
( gg7| ). Hence, we neglect the other just mentioned ad- 
ditive renormalizations. For our setup with the fixed ex- 



ternal field B, (3.47) implies for the free energy that 

(3.48) 



3~ rcn(^) * *F n 



(l),(2)_g2 



According to definition ( p. 73 ) this implies 

x (l).(2) (T , Wjff) ^(lH2) 

= X (1) ' (2) (r, w, u, M ) + x {lh(2) (r, w, u, M )div(3.49) 

for the susceptibilities. 

Now we are in good shape to analyze the scaling be- 
havior of the susceptibilities. In order to reduce the use 
of indices and to keep the arguments as plain as possi- 
ble, we carry out the following steps at the instance of 



X 



(i) 



At the end it will be straightforward to adapt our 
arguments to x^ ■ 

Just like the bare order parameter correlation func- 
tions, x^ 1 ) has to be independent of the inverse length 
scale (J,, i.e., it satisfies the identity 

d 

— vW (f,w,g) 



0. 



(3.50) 



This identity is now taken as the origin of an RGE for 
X^ ■ Expressing the bare quantities through their renor- 
malized counterparts one arrives initially at 



— + 13— 



d , d 

TK— + W( — 
OT OW 



(r, w, u, n) + x (1) (t, w, u, M)div} 



= 0. 



(3.51) 

Here it is important to realize that all the individual 
terms appearing in the RGE have to be free of e-poles. 
All terms associated with e-poles must cancel order by 
order in the loop expansion. T aking into account the 
form of the Laurent expansion (3.43), the form of the 

Wilson function (3 as given in ( 3.9b ) and that Xdw ~ L i ~ £ 
we obtain the RGE 



d_ 

dpi 
Tt 2 



■0 



d_ 

du 



8 



d_ 

dw 



X (1) (t,w,u,h) 



Xi{u)pL" 



(3.52) 
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where X\(u)fi e = (1 + ud u )Xi(u). 

bmce the RGE ( HH) is inhomogeneous, its solution is 
of the form 



X (1) (t, w, u, /x) = (t, w, u, m) + Xp 1 ' (t, w, u, /x) 



(1) 



(3.53) 



where xj, is the general solution of the corresponding 

homogeneous equation and Xp^ is a particular solution 
of the inhomogcncous equation. At the fixed point u* 
the method of characteristics gives for the homogeneous 
solution 

X^ ) (r,w,u,fi) = x ( k ) (l K *T,t?w,u;tn) , (3.54) 

where k* — n(u*) and C* = C( u *)- This solution has to 
be complemented by a dimensional analysis, 



Xh^ (t, w, u, fi) = n d 4 Xh 1 2 t,H 2 w,u,l) .(3.55) 



From Eqs. (3.54) and (3.55) we deduce that 



X^ (r,w,u, f x)=£ d - 4 x^ (r^r.r^tfl.ti*^) ■ 

(3.56) 

Exploiting our freedom to choose the flow parameter, we 
set I = \t\ v . This choice yields 

(t,w,u,h) = |r|( d - 4 >/ (1) (w\t\-*) (3.57) 



with /^^ being a scaling function. Due to Eq. (3.41) 
we know that ~ x 



Thus, the homogeneous 



solution may be written as 



t—2u 



(3.58) 



where t = (d — 2)v + <fi is the conductivity exponent of 
the RRN. 

A particular solution is readily found by making an 
ansatz that is as similar to the inhomogeneity as possible. 
We obtain 



xi 1} (t,w,u,h) 



x lfi - 



Tt 2 



w 2- e-2/u + <j>/v 



, (3.59) 



where X 1 = X\(u*). 

Combining Xh > Xp an( i the terms varying as a power 
of A gives us the scaling behavior of X ; vrz -i 



X 



(i) _ 



T 



A 4 + A 2 |r| 1 



T\t\ 



(3.60) 



where we have replaced nonuniversal coefficients, that 
might obscure the essential structure, by unity. Finally, 
we recall that our coarse grained parameter w is propor- 
tional to the temperature. Ignoring the proportionality 
constant we thus get 



X 



(i) 



A 4 + A 2 M 



\t-2v 



(3.61) 



The first terms correspond to the beginning of the small 
t expansion of the regular part of x ■ The last term 
characterizes the leading behavior of the singular part of 
X*- 1 ). Note that, over all, the leading small r behavior 
of x*- 1 ^ is determined by its regular part. In passing we 
mention that the coefficients of the |t| 2 term and the 
singular term form a universal ampl itude ratio that may 
be calculated with help of Eq. ( 3.59 ) 

An analysis similar to the one presented in the preced- 
ing paragraphs can be applied to x^ ■ Since only slight 
modifications are required, we simply state the result 



X 



(2) 

1 

T 

It 



A 



i 6 



A' 

2t-dv 



T 



A 2 M 



A 4 |t|+A 2 |t| 2 



[1 + T|r|-*] 



(3.62) 



where we have once more replaced unimportant constants 
by unity. As above, the regular part gives the leading 
small |r I behavior. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

In summary, we have studied the scaling behavior dia- 
magnetic susceptibility x^ and the mean square fluctu- 
ations of the total magnetic moment x^ 2 "* of large clusters 
as the Meissner (transition I) and the SG phase (transi- 
tion II) are approached as p — > p c at low T. Our main 
results are summarized by the formulae 



/(!) 



A roe 



\t — 2v 



= 



for transition I and 



A 1 ) 



X 



(2) 



0. 

y (2) 
A roe 



(4.1) 
(4.2) 

(4.3) 



\2t—dv 



T 



[1 + T\t\-+] (4.4) 



for transition II. We remind the reader that r measures 
the distance from the respective transition, r ~ p c — p. v 
is the percolation correlation length exponent and t and <p 
are the conductivity and the resistance exponent, respec- 
tively, of the RRN. Xreg and Xreg summarize the regular 
parts of the susceptibilities. These regular parts are very 
important and must not be neglected. In fact, they deter- 
mine the leading small r behavior of x^ and x'' 2 -' ■ As far 
as the leading small T behavior is concerned, our results 
capture anticipated features of the susceptibilities. Typ- 
ical for a diamagnetic susceptibility, x approaches at 
transition I for T — » a finite constant, x^ ° n the other 
hand diverges at transition II as T^ 1 . That is because 
X*^ represents a paramagnetic rather than a diamagnetic 
susceptibility due to randomly frozen magnetic momenta 
in the spin glass phase. 
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We point out that our results hold to arbitrary order 
in e-expansion. Using the powerful methods of renormal- 
ized field theory, we were able to explored general struc- 
tural properties of the Feynman diagrams contributing 
to the susceptibilities. This allowed us to determine the 
scaling behavior of x^ an d X to arbitrary order in 
perturbation theory. 

Our work reveals that the results by JLW are not en- 
tirely correct. On one hand, JLW overlooked the addi- 
tive character of the renormalization of the susceptibil- 
ities that leads to regular contributions. On the other 
hand, JLW did not realize that the coupling constant u>2 
associated with T 2 can enter the susceptibilities only in 
form of an irrelevant combination with w. This is true 
though and consequently exponents associated with T~ 2 
do not enter into the leading singular behavior of the 
susceptibilities. 

Closing, we would like to mention interesting issues for 
future work. For example one could extend our work by 
investigating the diamagnetism beyond transition I and 
II, i.e., for the other transitions featured in the phase 
diagram of the RJN. One might look at the normal to 
superconducting transition as p — > p c away from the im- 
mediate vicinity of T = [cf. the critical surface (1) in 
Fig. |J. Using our Gaussian propagator, this should be 
a feasible task. Another interesting subject is the role of 
vortex excitations, that we believe to be important in two 
dimensions. To address this question one has to develop 
an approach that avoids a linearization of the network 
equations. A potential strategy jll] is to devise a Villain 
type model @ for the RJN. 
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APPENDIX A: THE PROPAGATOR 

In this Appendix we derive the Gaussian propagator 
as given in Eq. (3.1). For notational simplicity we set 
T = 1. 

The task is to solve the differential equation 



r(A) - (V-iA- A 



G(x,x',A) = <S(x-x') .(Al) 



Here we use t he sho rthand notation r(A) = r+w\ 2 intro- 
duced in Sec. Ill A . Choosing for convenience the Landau 
gauge A(x) = B(0, x, 0j_) , Eq. ( [All ) takes on the form 

r(A) - W 2 + 2iuj(\)xd y +uj(X) 2 x 2 ] G(x,x',A) 
= £(x-x'). (A2) 



Inspired by Lawrie we rewrite the propagator as 

>(A) 



G(x, x', A) = cxp 
x G(x, x', A) , 



-(x + x')(y - y') 



(A3) 



where w(A) is the cyclotron frequency given in Eq. ( [2.79 ) 
Inserting the propagator ([A3]) into Eq. dA2) we get 



-(A) - V 2 + iw(A) [{x - x 1 ) d v -(y- y') d x ] 



;(X) S 



<5(x- 



(x - x'f + (y - y')' 



G(x,x', A) 



(A4) 



From Eq. flA4| ) we deduce that G(x, x', A) must be a func- 
tion of the difference of the coordinates x and x', i.e., 
G(x, x',A) = G(x — x',A). Hence, it is convenient to 
switch to momentum space via the Fourier transforma- 
tion 

G(x,A)= J G(k,A)exp(ik-x) (A5) 

with k = (p, q, kj_). Because the system is rotation- 
ally invariant in the x-y-plane as well as in the hyper- 
plane perpendicular to the x-y-plane, we anticipate the 
following form in momentum space: G(k, A) = G(jp 2 + 
q^k^A). Thus, Fourier transformation of Eq. (A4) 
leads to 



"(A) 



k 2 - 



[d 2 P + d; 



G(k, A) = 1 . (A6) 



For vanishing cyclotron frequency, this reduces to the 
well known equation for the propagator of the RRN. In 
Schwinger representation, the RRN propagator reads 



G(k, A) = J exp -s (V (A) + k 



(A7) 



However, our interest is not limited to the particular case 
ui(\) = 0. Hence, we generalize the solution (A7) by 



making the ansatz, 
G(k,A) - 



dsf(s) 



x exp 



-s{r(X)+kl)-g( S )(p 2 + q 2 )] (A8) 



with /(s) and g(s) being unknown functions of the 
Schwinger parameter s. We demand that 



r(A)+k 2 - 
{/(s)exp 



G(k,A) 



A 9 

as — 

OS 



(r(X)+k 2 ± )- 9 (s)( P 2 + q 2 )}} 



(A9) 
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For Eq. (A9) to be satisfied, the unknown functions f(s) 
and g(s) have to satisfy the differential equations 



<?'(*) =1- 



f(s) 



w(A)g(s) , 
= ^(A) 2 ff(s) . 



(A10) 
(All) 



along with the boundary conditions /(0) = 1, g(0) = 0, 
and g(s) > 0. And, of course, the two functions must 
yield f(s) -> 1 and g(s) -> for w(A) -> 0. We find 



tanh 



9(s) = ' 
/(«) = 



(u,(A>) 



w(A) 
1 



cosh ( w(A)s J 



(A12) 
(A13) 



Inserting these results into our ansatz (A8) we obtain 
G(k, A) as given in Eq. (fTtj). 
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